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Abstract 

Neutrino oscillations are only observable when the neutrino production, propaga- 
tion and detection coherence conditions are satisfied. In this paper we consider in detail 
neutrino production coherence, taking tt — >■ [iv decay as an example. We compare the 
oscillation probabilities obtained in two different ways: (1) coherent summation of the 
amplitudes of neutrino production at different points along the trajectory of the parent 
pion; (2) averaging of the standard oscillation probability over the neutrino production 
coordinate in the source. We demonstrate that the results of these two different ap- 
proaches exactly coincide, provided that the parent pion is considered as pointlike and 
the detection process is perfectly localized. In this case the standard averaging of the 
oscillation probability over the finite spatial extensions of the neutrino source (and 
detector) properly takes possible decoherence effects into account. We analyze the 
reason for this equivalence of the two approaches and demonstrate that for pion wave 
packets of finite width a xn the equivalence is broken. The leading order correction 
to the oscillation probability due to a xn ^ is shown to be ~ [ v g/( v g ~ v n )]cr X7T /l osc , 
where v g and v n are the group velocities of the neutrino and pion wave packets, and 
l osc is the neutrino oscillation length. 
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1 Introduction 

It is well known that neutrino oscillations are only observable when the conditions of coher- 
ent production, propagation and detection of different neutrino mass eigenstates are satis- 
fied. The production and detection coherence conditions ensure that the intrinsic quantum- 
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mechanical energy uncertainties at neutrino production and detection are large compared 
to the energy difference AEjk of different neutrino mass eigenstates: 



AE 3k ~ « a E , (1) 

where <je = niin{cx^ rod , o"^, et }. If this condition is violated, the neutrino production or 
detection process will be able to discriminate between different neutrino mass eigenstates, 
thus destroying neutrino oscillations. Indeed, neutrino mass eigenstates do not oscillate in 
vacuum, and the oscillations are only possible because neutrinos are emitted and detected as 
flavour states - coherent linear superpositions of different mass eigenstates. This coherence 
is destroyed if the condition (TTJ is violated. Production/detection coherence is known to be 
related to the localization of the corresponding neutrino emission and absorption processes, 
as it is the localization of these processes that determines the energy uncertainty <je [I]. 

The propagation decoherence can take place if neutrinos propagate very long distances. 
It is related to the fact that the wave packets describing the different neutrino mass eigen- 
states that compose a flavour state propagate with different group velocities, and therefore 
after a long enough time they cease to overlap and separate to such an extent that the am- 
plitudes of their interaction with the detector particles cannot interfere. For ultrarelativistic 
neutrinos the propagation coherence condition can be written as 

~2]0- L < ® xv - v 9 /°e , (2) 

where L is the baseline, v g is the average group velocity of the wave packets of different 
neutrino mass eigenstates and a xv is their common effective spatial width. Note that both 
the production/detection and propagation coherence conditions, eqs. (CEJ) and (j2J), put upper 
limits on the mass squared differences Am| fc of different neutrino mass eigenstates. 

In practice, the propagation coherence condition is very well satisfied in all cases except 
for neutrinos of astrophysical or cosmological origin, such as solar, supernova or relic neu- 
trinos. As to the production/detection coherence condition, it is usually tacitly assumed to 
be always satisfied due to the extreme smallness of the neutrino masses (and therefore of 
their mass squared differences). However, there may exist situations in which this is not the 
case. In particular, in view of possible existence of sterile neutrinos with masses in the eV, 
keV or even MeV range, the question of production/detection coherence should be analyzed 
especially carefully when the corresponding mass squared differences are involved. Note 
that there are some hints for such sterile neutrinos coming from short-baseline accelerator 
experiments, the reactor neutrino anomaly, gallium radiative source experiments, r-process 
supernova nucleosynthesis, pulsar kicks, warm dark matter and leptogenesis scenarios [2]. 

The issue of coherence is of utmost importance for neutrino oscillation experiments. If 
the coherence conditions are strongly violated, the probabilities of flavour transitions will 
correspond to averaged out oscillations, i.e. they will have neither L nor E dependence. This 
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means that the two most important signatures of neutrino oscillations - the dependence on 
the distance between the neutrino source and detector and the distortion of the energy 
spectrum of the detected signal - will be absent. In particular, the two-detector setups of 
neutrino oscillation experiments will be completely useless. In such situations one would 
have to rely entirely on the overall normalization of the neutrino flux, which is usually known 
with insufficient accuracy. 

The question of neutrino production coherence has been recently addressed in ref. [3], 
where it was considered in the framework of incoherent summation of probabilities of neu- 
trino production at different points inside the neutrino source. In addition, a simplified study 
of the coherent amplitude summation was performed in that paper. It was demonstrated 
that under certain conditions the two summation procedures lead to identical results. 

In this paper we address this question within a more rigorous and consistent approach. 
We study in detail neutrino production coherence, taking tc — > fiu decay as an example. 
With minimal modifications, our analysis will also be applicable to neutrino detection. We 
find the oscillations probabilities in two different approaches: 

1. An approach based on the quantum-mechanical wave packet formalism. We first cal- 
culate the transition amplitude by summing the amplitudes corresponding to neutrino 
production in pion decay at different points along the trajectory of the parent pion. 
We then calculate the transition probability and study its coherence properties. 

2. We assume that for each individual neutrino production event the oscillation probabil- 
ity is fully coherent but depends on the exact position of the neutrino production point, 
i.e. is described by the standard oscillation formula. The effective oscillation proba- 
bility is then found by integrating (averaging) this standard coordinate-dependent 
probability along the neutrino source, with the proper exponential factor describing 
the pion decay included [3]. 

It should be stressed that, strictly speaking, the Erst approach, based on a quantum me- 
chanical amplitude summation, should always be used. The second (probability) summation 
procedure, which is classical in nature, is however much simpler and is usually employed in 
the analyses of experiments. It is therefore one of the main goals of the present paper to 
study if and when the use of the probability summation approach is indeed justified. 

We compare the results of the two approaches and demonstrate that they exactly coin- 
cide when the parent pions are considered as pointlike particles and the detection process is 
perfectly localized in space and time. One can therefore conclude that for pointlike parent 
particles and well-localized detection the standard averaging of the oscillation probabilities 
over the finite spatial extensions of the neutrino source (and detector) properly takes possi- 
ble decoherence effects into account. We analyze the reason for this equivalence of the two 
approaches and demonstrate that for finite-size pion wave packets the equivalence is broken. 
We show that for small widths of the pion wave packets <y xn the oscillation probabilities get 
small oscillatory corrections which are linear in a xn . For different shapes of the pion wave 
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packets these corrections have the same form and differ only by numerical coefficients. At 
the same time, large lead to production decoherence and thus to suppression of 

the oscillations. We also consider the production coherence in the case when the charged 
leptons accompanying neutrino production are detected, leading to neutrino tagging, as well 
as in the case when the interactions of pions in the bunch between themselves or with other 
particles which may be present in the neutrino source are taken into account. 

A note on terminology. In this paper we use the word 'coherence' in two different, 
though related, senses. By production coherence we mean that different neutrino mass 
eigenstates are produced coherently, so that the emitted neutrino is a flavour eigenstate. By 
coherent summation we mean the summation of the amplitudes of neutrino production at 
different points along the trajectory of the parent pion. As we shall see, these two coherences 
are in fact acting in opposite directions: if the production is incoherent, then the coherent 
summation (i.e. approach 1 discussed above) is mandatory, whereas for coherent neutrino 
production the incoherent probability summation (approach 2) is justified. The common 
feature of the two coherences is that both require summation of certain amplitudes. Indeed, 
in the case of neutrino production coherence, these are the amplitudes of emission of different 
neutrino mass eigenstates, whereas in the case of coherent summation approach, these are 
the amplitudes of neutrino emission from different space-time points. 

The paper is organized as follows. In section [2] we calculate the oscillation probability 
for neutrinos produced in decays of pointlike pions in a decay tunnel. The calculations 
are performed in the quantum-mechanical wave packet approach with summation of the 
amplitudes of neutrino production at different points along the tunnel. In section [3] we 
recapitulate how the same problem is solved at the level of summation of probabilities 
rather than amplitudes. In section H] we generalize the results of the amplitude summation 
approach to the case of non-zero spatial width of the pion wave packets. We then consider 
the particular cases of Gaussian and box-type pion wave packets. In section [5] we apply the 
probability summation approach to the case of protons incident on a finite-thickness target. 
We also compare the obtained results with those of section |H In section [H] we develop an 
alternative approach to the quantum mechanical amplitude summation calculation, which 
is valid for arbitrary shapes of the pion wave packets. In section [7] we consider effects of 
possible detection of the charged lepton accompanying the neutrino production on neutrino 
production coherence. We also briefly discuss here the effects of the interaction of pions in 
the bunch between themselves or with other particles which may be present in the neutrino 
source. Section [S] is devoted to implications of our analysis for various neutrino experiments. 
We summarize and discuss our results in section [91 
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2 Oscillation probabilities in the wave packet approach. 
Pointlike parent particle approximation 



In the quantum-mechanical wave packet approach, the oscillation amplitude A a p(L,t) is 
obtained by projecting the evolved neutrino state, which was initially produced as the flavour 
eigenstate z/ Q , onto the detected neutrino flavour eigenstate vp. The oscillation probability 
P a p{V) is then found by integrating the squared modulus of A a ^{L,t) over time (see, e.g., 

ID)- 

Let us first consider oscillations of neutrinos produced in decays of quasi-free parent 
particles. By this we mean that the decaying particles may be confined to a finite-size 
source or a decay tunnel, but their interactions with each other or with other particles 
which may be present within the source can be neglected (we will relax these assumptions 
in section [7]). As an example, we take the tt — » decayQ We will treat this problem 
in a 1-dimensional approach, i.e. assuming that the neutrino momentum p is parallel to 
the baseline vector L connecting the pion production and neutrino detection points. Such 
an approximation is well justified when the transverse (i.e. orthogonal to L) sizes of the 
neutrino source and detector are small compared to L. 

To find the wave function of the produced neutrino state, we will need the wave functions 
of the parent pion and the muon which participate in the production process. 



2.1 The pion and muon wave functions 

We will be assuming that pions are produced by a beam of protons incident on a solid- 
state target and then decay inside a decay tunnel, with the total length of the decay tunnel 
being l p . Such a setup corresponds to accelerator experiments. The nuclei in the target are 
well localized, with the uncertainty of their position <j x n being of the order of inter-atomic 
distances, <j x n ~ 10~ 8 cm. The spatial width of the wave packets of the incident protons a xp 
depends on the conditions of their production; it cannot exceed the mean distance between 
the protons in the bunch, and e.g. for the Fermilab NuMI source it can be estimated as 
a xp < 1CT 4 cm. The spatial size of the pion production region a x p can be defined as (6J [7J 

111 

+ T-T- (3) 



2 

a xP a xN a xp 

We can also define the effective velocity of the pion production region as the sum of the 
group velocities v p and vn of the incident proton and the target nucleus, weighted with the 
inverse squared widths of the corresponding wave packets: 

(4) 

V a xp a xN J 



1 A similar (though somewhat different and less detailed) analysis was performed in [5]. 
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The coordinate-space width of the pion wave packets can then be found from the relation 

a xn « [a 2 xP + (v P - VvY/aU 1 ' 2 , (5) 

where v % and <r en are, respectively, the group velocity of the pion wave packet and the energy 
uncertainty of the produced pion state. The latter is approximately equal to the inverse of 
the overlap time of the proton and nucleon wave packets at pion production [7]: 

V P -V N Vp_ 

Eq. (E]) thus has a simple physical meaning: the first term in the square brackets is the 
contribution of the finite spatial size of the pion production region to the width of the pion 
wave packet, whereas the second term is related to the fact that the pion production takes 
finite time. 

As follows from ([3]), the coordinate uncertainty of the pion production point a x p is 
dominated by the size of the shortest wave packet of the participating particles, which in 
our case is a X N- cr x p ~ cf x n- We shall also assume that the target nuclei are at rest, vn = 0. 
Eq. (j3J then gives vp ~ v p (a X N / '<J xp ) 2 <C v p . From eqs. fl5]) and (EJ) we find 

Kn + «v p ) 2 a 2 xp ] 1/2 ■ (7) 

Thus, a xn < a xp < 10~ 4 cm, i.e. the spatial width of the pion wave packet is much smaller 
than all the lengths of interest in the problem - the length of the decay tunnel l p , the baseline 
L and the oscillation length Z osc = A7cp/Am 2 , where p is the neutrino momentum. Therefore 
to a very good approximation one can consider the pions as pointlike particles. 

Let us recall that the coordinate-space wave packet describing a moving free particle 
can be written, in the approximation where the spreading of the wave packet is neglected, 

as 

where p is the peak momentum of the wave packet, v g = (dE(p)/dp)\ P0 is its group velocity 
and g(x — v g t) is its shape factor (envelope function) (see, e.g., [9]). The shape factor is the 
Fourier transform of the momentum distribution amplitude (i.e. of the momentum-space 
wave function) of the pion; it quickly decreases when \x — v g t\ becomes large compared to 
the spatial width of the wave packet o x . Since g depends on x and t only through the 
combination x — v g t, for stable particles eq. §E§ describes a wave packet propagating with 
the group velocity v g without changing its shape. For unstable particles the energy E(p ) 
in eq. ([8]) should be replaced according to E(p ) — > E(p ) — iT/2, where T is the particle's 
decay width in the laboratory frame. 

As was pointed out above, in the problem under consideration pions can to a very good 
approximation be regarded as pointlike particles, i.e. the shape factor g n (x — v^t) of their 
wave packets can be taken to be a 5-function. The pion wave function then takes the form 

fa(x, t) = a e «J*-U5-«)t-rt/2 s( x - v n t) box(x; l p , 0) , (9) 
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where C w is a normalization constant, Q is the mean momentum of the pion state, E n (Q) = 
(Q 2 + m 2 ) 1 ^ 2 , v n is the group velocity of the pion wave packet, and the function box(;c; A, B) 
is defined as 

box(*;AB) = {j; l^i/' • (10) 

The box function in ([9]) takes into account that pions are produced at the beginning of the 
decay tunnel and that undecayed pions are absorbed by the wall at the end of the tunnel. 
The pion is assumed to have been produced at the time t = 0, and the factor e~ vt l 2 where 
T is the pion decay rate in the laboratory frame takes into account the exponential decay 
of the pion's wave function. Obviously, the production at t = x = is an approximation, 
as the pion production process has finite space-time extension. Strictly speaking, eq. is 
only valid outside the pion production region, whose size we neglect here. 

In general, when the spatial width of the pion wave packet a X7T is considered to be finite, 
one can also employ the pion wave packet of the type with S(x — v n t) replaced by the 
corresponding finite- width shape factor g w (x — v n t), provided that the effect of the pion 
production process itself on neutrino oscillations can be neglected. This approximation is 
valid when the spatial width of the pion wave packet is negligibly small compared to the 
neutrino oscillation length / osc and, in addition, the pion decay effects during its formation 
time are negligible: 

Ta X7r /v n < 1 . (11) 
We shall consider the case of finite a X7T in sections H] and [6j 

Throughout most of this paper, we shall be assuming that the muon produced alongside 
the neutrino in the pion decay is undetected and that its possible interaction with the 
environment can be ignored. In this case the muon is completely delocalized (a Xfl — > oo) 
and therefore can be described by a plane wave [8]. The effects of possible muon detection 
or interaction with medium will be considered in section [7J 



2.2 Calculation of the neutrino wave packet 

To find the neutrino wave packet, we first calculate the amplitude A(n — > fii/j) of the pion 
decay with the production of a mass-eigenstate neutrino Vj of mass rrij. We will be assuming 
that the pion is described by the wave function ([§]), while the muon and the neutrino are 
described by plane waves of momenta K and p, respectively. This will give us the probability 
amplitude that Vj and the muon are produced with the momenta p and K. The standard 
calculation gives 

/oo /*oo 
dt / dx e i[EAP)+E,(K)]t-i( P+ K) X t) _ (12) 

oo J — oo 

Here Ej(p) = (p 2 + m 2 ) 1 ^ 2 , E^(K) = [K 2 + m 2 ) 1 / 2 , and Mp is the coordinate-independent 
part of the pion decay amplitude. In obtaining ffTiZl) it was taken into account that the pion 
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decay amplitude Mp(q,k) is a smooth function of the pion and muon momenta q and k 
and therefore it can be replaced by its value Mp = Mp(Q, K) taken at the mean momenta 
(for the muon, k = K). The leptonic mixing parameter is not included in the Vj 
production amplitude (1121) since it will be explicitly taken into account in the definition of 
the oscillation probability. Though the spatial integration in ( IT2l) is formally performed in 
the infinite limits, the box function in the expression (Q for ip n (x,t) implies that in reality 
it only extends over the interval [0, l p \. From (Q it also follows that the integral over time 
in f|T2|) receives non-zero contribution only from the region t > 0. 

For a fixed K the quantity f^(p,K) gives the amplitude of the neutrino momentum 
distribution, i.e. the momentum-space wave packet fj(p) of the produced Substituting 
(|9]) into ( 1T2|) and performing the integrations, we obtain 

1 _ e i[Ej{p)-Ep-v 7r ( P -P)+ir/2]l p /v^ 

f > {p) = Cl E 3 {p)-E P -v^p-P) + rT/2 ■ (13) 

Here Cj is a constant, and the quantities P and Ep are defined as 

P = q_ K: E P = E V (Q) - E^K) , (14) 

i.e. in the plane-wave limit they would be the momentum and energy of the emitted neu- 
trino. The momentum distribution amplitude f?(p) contains the usual Lorentzian energy 
distribution factor corresponding to the decay of an unstable parent state. The second term 
in the numerator of (TT5|) (the exponential phase factor) reflects the fact that the parent pion 
is not completely free but is confined to a tunnel of length l p ; it would be absent in the case 
of decay of free pions, which can be considered as the limit l p — > oo. 

The oscillation amplitude A a /3(L,t) can now be obtained as a projection of the evolved 
neutrino state onto the detected one directly in the momentum space. However, the 
coordinate-space approach is more illuminating, and therefore we present it first. Momentum- 
space calculations will be employed in section 0] for the case a xn ^ 0. 



2.3 Neutrino wave packet in the coordinate space 

The neutrino wave function in the coordinate space ipj(x,t) can be obtained by Fourier- 
transforming the momentum-space neutrino wave packet f|T3l) : 

^(x, t) = / — ff (p) e -*%(p)*-H»* . (15) 
J J 2tx J 

Substituting f fT3|) into f lT5|) . expanding Ej{p) near the point p — P up to terms linear in 
p — P and performing the integration, we find 

^(x, t) = const. e-i^P^+iP^L-WJ^)^-^ box (^ _ x . V -l^^l p) o) } . (16) 
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Figure 1: Schematic representation of the coordinate dependence of the neutrino wave packet 
( TT5|) (arbitrary units). 



Here the following notation has been used: 



Ps = P 



E P ~E 3 {P) 



where 



Vj - V-n 



V 3 



E s {Pi) - Ej{P) + ,•(/> - P) = /;,(/') + v. 



_dEj(p), P 
~~ dp 'p= p ~ Ej(P) 

Ep-E 3 (P) 



(17) 



Vj - v„ 



Note that for finite l p the momentum distribution \fj{p)\ 2 is an oscillating function of p, 
and p = Pj corresponds to the peak of the envelope of this function. 

The wave packet in eq. (TT6T) has the general form (jSJ) with the shape-factor function 



g(x- Vjt) 



2(vj 



boxrvjt 



Vj - V 



-l P , o) 



(19) 



The presence of the box function here means that the neutrino wave packet has sharp front 
and rear edges 0. The box function enforces that the front of the neutrino wave packet arrives 
at the point x at the time t\ = x/vj and leaves this point at t 2 = x/vj + (l/v n — l/vj)l p u 
The distance between the edges of the wave packet is thus Vj(l/v n — l/vj)l p . The shape 
factor reaches its maximum at the front of the wave packet and decreases exponentially 
towards its end. This gives an interesting example of an asymmetric coordinate-space wave 
packet (see fig. [TJ. 

Having found the expression for the neutrino wave function, we can now calculate the 
expectation value x of the neutrino coordinate as well as the coordinate dispersion, which 
gives the spatial width of the neutrino wave packet a xv . In the limit Tl p /v n ^> 1, when the 



2 Note that this is related to our assumption of pointlike pions. For pion wave packets of finite spatial 
size (and e.g. Gaussian form) the neutrino wave packets would not have sharp borders. 

3 This has a simple interpretation: the emission of neutrino wave packet by the pion starts at t — and 
abruptly ends when the pion reaches the end of the decay tunnel, i.e. at t = Ip/v^. Therefore t\ — x/vj and 



h = Ip/Vir + (X— l p )/Vj = x/Vj + - 1/Vj)l 



V 
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length of the decay tunnel l p is large compared to the pion decay length /decay = v n /T, we 
find 

Vj - V n Vj - v n 

x « Vjt - - J — — , o xv « — . (20) 



Note that in this limit the box function in eqs. ( I16p and ( II 9p has to be replaced by 9(vjt — x 
where 6{x) is the Heaviside step function. 

In the opposite limit, Ylp/v^ C 1, we find 

Vj — V n 1 Vj — v 



-~ Vj - v n , ^ J- <Jj - Vtt j / 9 -,s 



Notice that in this limit only a small fraction of pions decays before being absorbed by the 
wall at the end of the decay tunnel. 

Although, as was mentioned above, the distance between the sharp borders of the neu- 
trino wave packet is [(vj — v n )/v n )lp, the effective spatial width a xv of the wave packet is not 
necessarily determined by this quantity: this is only the case when the exponential decrease 
of the envelope function is relatively slow, Tl p /v n <C 1 (see (jUJ). In the opposite case, 
Tlp/v-x ^> 1, which corresponds to the decay of unconfined free pions, it is this exponential 
decrease that dominates a xu , as described by eq. ( 12"U|) . It is interesting that the neutrino 
wave packets are of macroscopic size in the case we consider. 

Next, we have to find the amplitude Aj(L, t) which is the contribution of the jth neutrino 
mass eigenstate to the transition amplitude ^4 Q/ g(L,t). In the coordinate-space calculation 
it is simply given by the projection of i/jj(x,t) onto the detected neutrino state ^P(x): 



Aj(L, t) = j dx rl>?*(x) ip°(x, t) . (22) 

The detected state is in general a wave packet centered on the point x = L. In what follows 
we will be assuming that the detection process is well localized both in space and time, so 
that the wave function of the detected neutrino state can be represented in the configuration 
space by 8{x — L). This requires some explanations. The momentum-space wave function 
of the detected neutrino state fP(p) is given by a formula similar to that in eq. (fl2|) . with 
the integration over the detection coordinate (£', x'). We assume that the detection process 
is perfectly localized in space and time, so that the integrand of this formula contains the 
factor 6(x'—X£))S(t , —t£)). As a result, fP(p) is actually momentum independent, and for the 
coordinate-space wave function of the detected neutrino state centered on the point x = L, 
which is a time-independent Fourier transform of ff{p) [U [8], we find 4>f{x) = S(x — L). 
From (T22]) we then obtain 

A J (L,t)=rf(L,t). (23) 
The oscillation probability P a p{L) is given, up to a normalization factor, by (see, e.g., jl]) 

P aP {L) = J2 KjUpjU ak U; k I jk (L) , (24) 
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where 



From eq. fl23l) we then find 



I jk (L) = / dtAj(L,t)Al(L,t) 



(25) 



Ijk(L) 



dt^(L,t)i; s k *(L,t). 



(26) 



Substituting here the expressions for ipj(x) and ip^( x ) fr° m ( tlSD and performing the inte- 
gration, we obtain 



2P 



e 



Am 2 , Am 2 , 

£ _ p -n P /^ e -i^(L-« P ) 



(27) 



Here we have discarded terms ~ (Amjj.) 2 and ~ Am 2 k T/P and neglected the difference 
between v j and v k whenever they are multiplied by Y or by Am The constant 

C\ can be found by imposing the unitarity constraint 'Yl l aP a p{L) = 1 on the oscillation 
probability ff2^|) FI which leads to the normalization condition Ijj(L) = 1 |4J. This yields 
C% — (1 — e~ rip / V7 ')~ 1 , so that we finally obtain 



Ijk(L) 



iV 



'l_ e -rw«-) H 



+ ^^ 



e 2P 



Am 2 j, Am 2 , 



(28) 



By making use of the normalization condition Ijj(L) = 1 one can demonstrate that the 
quantity Ijk{L) (and therefore the oscillation probability (pM)) ) is Lorentz invariant [I]. 



2.4 The oscillation probability 

With the expression for Ijk{L) at hand, we can now calculate the oscillation probabilities 
from eq. (El]). Consider a 2-flavour case when the oscillations are governed by just one mass 
squared difference and one mixing angle. As an example, short-baseline oscillations in the 
3+1 scheme with one sterile neutrino are essentially reduced to effective 2-flavour ones with 
Am 2 = Am^j ^ 0, while the much smaller mass squared differences Am^ and Am^ can be 
neglected. In this case the survival probabilities P aa are described by the effective mixing 
parameters 

wn.6 = 8=\U oA \, cosfl = c= (1 - \U a4 \ 2 ) 1/2 . (29) 

4 Note that the latter approximation implies that we neglect the effect of decoherence due to the wave 
packet separation (propagation decoherence). This is justified for L <C l co h, where Z C oh — (v g /\vj — Vk\)o~ X v, 
which is the case we are mainly interested in. The propagation decoherence effects can, however, be readily 
taken into account. 

5 This normalization prescription looks rather ad hoc in the quantum-mechanical wave packet approach 
but can be rigorously justified in the quantum field theoretic framework [5]. 



12 



In the case of transition probabilities P a p (a ^ (3), the mixing parameter can be chosen as 

sin 2 26 = A\U a4 \ 2 \U,3 4 \ 2 . (30) 

Consider, for example, the survival probability of muon neutrinos in the 2-flavour 
scheme. From eqs. ( l24l) and ( |28l) we find 

2c 2 s 2 1 

Pw = c 4 + s 4 + ^ 2 - i _ e _ rip/v ^ [cos + £ sin - e" r ^ [cos(0 - <j> p ) + £ sin(0 - p )]] . 
Here 

Am 2 Am 2 



(31) 



L ' 0P = ~n~E)~^P ' ( 32 ) 



2P ' r * 2P 
and we have defined the parameter 

Am 2 . , 

(33) 

which, along with Tl p /v n , characterizes possible decoherence effects at neutrino production 
(see the discussion in section l2~5l) . In the limit v n = 1 eq. ( }3~TT) coincides with result found 
in ref. [3]. 

We shall now analyze the probability ( 13T]) in the light of neutrino production coherence. 
2.5 Coherence violation at neutrino production 

Let us start with a qualitative analysis of neutrino production coherence. Recall that the 
production coherence condition requires that different neutrino mass eigenstates forming a 
flavor neutrino state be emitted coherently in the production process. This condition can 
be written as 

A£«a E , (34) 

where AE is the energy difference of different neutrino mass eigenstates and <je is the 
quantum-mechanical energy uncertainty inherent to the production process. As was pointed 
out in [I], for decays of non-relativistic parent particles contained in a box of finite size, ge is 
determined by the larger of the two quantities: the particle's decay width T and the inverse 
time between two subsequent collisions of the particle with the walls of the box. In our 
case, the "collision time" is the time l p /v n between the pion production and absorption of 
undecayed pions at the end of the decay tunnel. In general, for relativistic parent particles 
of velocity vp this non-relativistic energy uncertainty should be multiplied by the factor 
Vg/(v g — vp). This takes into account that the neutrino is emitted by a moving particle 
in the forward direction, so that the parent particle "chases" the produced neutrino wave 
packet [12l H] • We thus find that in our case 

v 

(je ~ max{r, iv/L} — . (35) 

Vg ~ V w 



13 



It is known that the energy uncertainty oe also determines the spatial width of the 
neutrino wave packet [I]: a xv ~ v g /<jE- From ( )35|) we therefore find that in the limiting 
case Tlp/v-x ^> 1 the width of the wave packet a xv ~ [v g — v n )/T, whereas for Ylp/v^ <C 1 
one finds o xv ~ [(v g — v^)/v^l p . These results are in full accord with our previously found 
values of a xv (see eqs. ( 1201) and ( 12T1) ). 

As follows from (TT 



AE=\E, (Pi )-E A P h )\^^- y (36) 

Taking into account eq. (15o^) . for the (de)coherence parameter AE/<je we then find 
AE Attt? A B 

~ Z p = p , (n p /v w <l); ~£ (n p /^>l). (37) 

0"£ 2_P erg; 

Let us now consider the probability (13T]) from the point of view of the production 
decoherence effects (see also [3]). As was mentioned above, these effects depend in general 
on two parameters, £ and Tl p /v n . Note that these quantities have different dependence on 
the pion lifetime T, and that there is a relation between them and the phase <f) p : 

£-^ = P - (38) 

The three quantities in this equation can also be expressed through the pion decay length 
^decay = v v /Y and the neutrino oscillation length Z osc = AnP/Am 2 : 

C p„ ^decay ^\ lp , „ lp , qq n 

4 = 2ir- , = , <p p = 2tt— . (39) 

'osc 'decay 'osc 

Consider first the limit Tl p /v n C 1. As follows from (1371) . the decoherence parameter 
in this case is <p p . Thus, one can expect strong decoherence effects for <p p ^> 1 and no 
decoherence in the opposite limit. 

The case Tl p /v n <C 1 is actually more easily studied by going to this limit in the 
expression for Ijk(L) and then calculating the oscillation probabilities rather than by directly 
expanding the expressions for oscillation probabilities. Taking the limit Tlp/v^ <C 1 in (1281) . 
we obtain 



Am % i 



Am 2 -, Am 2 

l_ 

2P 



J±L _ p -i^lH(L-l p ) 

2P~ l V 



(40) 



Note that this quantity satisfies the correct normalization condition Ijj(L) = 1 provided 
that Ato|j is understood as the limit Am^ — > 0. Substituting (1401) into (1241) . we find 

Pw = c 4 + s 4 + [ sin - sin(0 - p )] . (41) 
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In the limit P ^> 1 this expression yields P MM ~ c 4 + s 4 , which corresponds to averaged 
neutrino oscillations, whereas for p 1 it gives the standard oscillation probability 

pstand = c 4 + ^4 + ^ ^ ^ ^ 

as expected. Let us stress that the latter result does not depend on whether £ is small or 
large, though £ must still satisfy the relation in eq. (138]) . As follows from the same relation, 
in the case 4> p ^> 1 (strong decoherence) £ is automatically very large in the Tl p /v n <C 1 
regime. 

Consider now the limit Ylp/v^ ^> 1, when most pions decay before reaching the end of 
the decay tunnel. In this case, according to ( 1371) . decoherence effects should be governed by 
the parameter £. As expected, in the limit £ — > eq. fl3T|) yields the standard oscillation 
probability f|42l) . whereas for £ ^> 1 the oscillating terms in fl3T|) are suppressed and one 
obtains the averaged probability P MAt ~ c 4 + s 4 . 

Thus, we conclude that in the case of sufficiently small l p , when both Tl p /v w <C 1 
and (j) p -C 1, the production coherence condition is always satisfied (and the oscillation 
probability takes its standard form) irrespective of the value of £, which is irrelevant in 
this case. If Tl p /v n <C 1 but <p p > 1 (which implies l osc /2-n < / p <C /decay), the production 
coherence is moderately or strongly violated, depending on the value of (j) p . In the case 
TZp/tv ^> 1 the production coherence depends on the parameter £: it is satisfied for £ 1 
and strongly violated in the opposite case. 

Irrespectively of the value of Tlp/v^, strong production coherence violation always im- 
plies £ 3> 1, although large values of £ do not necessarily imply coherence violation. It is 
easy to see that the same applies to <p p . Therefore, in order to find out if the production 
coherence is violated one has to check the values of any two out of the three parameters 
rip/vn, £ and 4> p (the third one will then be given by eq. fl38|) ). It is convenient to choose 
these to be £ and (fi p . 

Indeed, from the above considerations it is easy to see that if ri p /v n is very small or 
very large, the production coherence is strongly violated in the case £ ^> 1, (p p ^> 1 and is 
satisfied in all other limiting cases (large £ and small <fr p , small £ and large <f) p , and small £, 
small <p p ). If Tl p /v n ~ 1, the parameters £ and P in the limiting cases are either both small 
or both large, and from eq. (|3T|) it again follows that the production coherence is strongly 
violated if £ 3> 1, 4> p S> 1- It is satisfied in the opposite case. 

In the intermediate cases one can expect moderate violation of the production coherence, 
which should lead to noticeable deviations of the oscillation probabilities from their standard 
(i.e. predicted under the coherence assumption) values. 
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3 Incoherent probability summation approach 



Let us now consider oscillations of neutrinos produced in pion decays inside a decay tunnel 
in a entirely different approach [3]. We will be assuming that each individual neutrino 
production event is completely coherent, so that the oscillations of the produced neutrino 
are described by the standard probability (I42p . At the same time, one has to take into 
account that the neutrino production can take place at any point along the decay tunnel, 
though the pion flux decreases exponentially with the distance from the pion production 
point. The oscillated neutrino flux at a detector at the distance L from the beginning of 
the tunnel is then given by 

F,(E, L) = F W (E, 0)r f P e~% P^ nd (E, L - x)dx , (43) 

Jo 

where F W (E, 0) is the initially produced pion flux. The unoscillated neutrino flux F® (i.e. 
the flux in the absence of the oscillations) is given by eq. (I43p with P^ nd (E, L — x) in 
the integrand replaced by 1. One can now define the effective oscillation probability as the 
ratio of the oscillated and unoscillated fluxes: 

Pf,(L,E) = F,(L,E)/F°(L,E). (44) 

The calculation is straightforward, and the result turns out to coincide exactly with the 
expression for P^(E,L) given by eq. (l3Tj) . Thus, the summation of the amplitudes of 
neutrino production at different points along the pion path performed in section [2] and the 
summation of the corresponding probabilities carried out in this section lead to the same 
oscillation probability. We will discuss the reason for this intriguing coincidence and the 
conditions under which it is broken in section [6] and in the Discussion section. 

4 Amplitude summation: The case of finite- width pion 
wave packets 

Let us now relax the assumption of pointlike pions that we employed in section [2j For the 
pion wave function we will use the expression similar to Q), but with 5(x — v n t) replaced 
by a general shape factor g^{x — v n t): 

^{x, t) = a e iQx - iE ^ Q)t - Tt/2 g n (x - v*t) box(x; l p , 0) . (45) 

We will be assuming that g n (x — v n t) is peaked at or near the zero of its argument and 
rapidly decreases when \x — v^t\ > a X7r , but otherwise will not specify the exact form of this 
function. Substituting fj4"5l) into f|T2|) . for the momentum-space wave packet of the jth mass 
eigenstate component of the produced neutrino we obtain 

ff{p) = M P C n f P dx jlEi(P)-E P ]%~i(P-P)*-^ - Ep) , (46) 

Jo 
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where g w (Ej(p) — Ep) is the Fourier transform of g w (x — v n t) e 2v ^ ^ hct 

/oo 
dt > S E ^ E ^' [g^-vj)e^\. (47) 
-oo 

Here the integration variable is t' = t — x/v^. Note that the integration has been extended 
to the negative semi-axis of t' here; this is justified only if for large negative t' the function 
g^—v^t') decreases sufficiently rapidly, and in any case faster than e~ r ^'^ 2 . This condition 
is satisfied for both box-type and Gaussian pion wave packets which we will consider below|§ 

Performing the integration over x in (1461) . we obtain 

1 _ AE 3 {v)-E P -vA.P-P)+i^ 

= Q £ , (p) _ £p _^_p ) + 1 r ^W " E ^ • (*) 

where Cj is a constant. This expression differs from the corresponding one in the case of 
pointlike pions ( IT5|) by the extra factor g n (Ej(p) — E P ). 

Let us now calculate the transition amplitude and the quantity Ijk(L). We will be again 
assuming that the detection process is well localized in space and time. The contribution 
of the jth neutrino mass eigenstates to the oscillation amplitude is then given, according to 
eqs. 023) and (EES]), by 

ML,t) = J ^//(p)e^-^W*. (49) 
Substituting this into eq. ( |25l) and performing the integration over t, we find 



W)= I ^^//(Pi)/f(P2) e i(w - p2)i 27r5[^( Pl )-^(P2)] 



Am" 



(2vr) 2 Jj 

a Jk L 1 f dp 

Here we have used 



e~^ L - I J 1 ffipi)ft (Pi + A™?*/ 2P) . (50) 



SlEjipx) - E k (p 2 )} ~ 5[^( Pl - p 2 ) + Am* fc /2£] , (51) 
which follows from Ej{p\) — E k (p 2 ) ~ (dE/dp)Ap + (dE/dm 2 )Am 2 . Substituting (T4H!) into 



"There is an ambiguity in the way the decay exponential is introduced in the pion wave function, which 
is related to the fact that for a X7T =/= one should take into account the finite duration <7 tjr of the pion 
production process. The decay exponential can be written as e~ r ( t ~~ t °^ 2 , where to is e.g. the time when the 
pion formation is completed, or the time when the pion production process starts. This ambiguity, however, 
plays no role if IV t7r ~ Tux^/v^ <C 1, which we assume throughout this paper (see eq. Note that 

the extra factor e rt °/ 2 can be absorbed into the normalization constant CV of the pion wave function and 
therefore does not affect the oscillation probabilities. 
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50]) . we obtain 



V 9 



2tt Ej(p) - E P - v n (p - P) + z| 

-,-riip_-, , / „ / Am 2 



^ _ e -i[E k ( P+ A m 2 j j2P)-E P -v^p+A m %/2P-P)-iL]^^ E ^ p+ Ag^j _ ^ x 
* ,_, / Am 2 , \ ,_, / Am 2 , _n . r ' (^^) 

E k (p+^) -E P -v w (p+^-P) -iE 

As usual, the constant Cj should be found from the condition Ijj(L) = 1. 

We shall now consider two special cases, Gaussian and box-type pion wave packets. 



4.1 Gaussian shape factor of the pion wave packet 



In this case 



which gives 



_ (x — v-nt) 

g v (x - v w t) = e ^ 



g{Ej(p) - Ep) = const, e < 2 



(53) 



(54) 



Let us calculate Ijk{L) in the limit Tl p /v n 3> 1, which corresponds to the situation when 
most pions decay before they reach the end of the decay tunnel. Substituting f }54"|) into (1521 
and expanding Ej(p) near the point p = P up to terms linear in p — P, we find 



Ijk{L) — Cj 



iter 

-*^IF- L ( 2<x 2 



r "g 

2(o 9 -%) U7r 



-iB. 



erfc 



+e 



20-2 "» r 



-iB. 



erfc 



r 



2(Vg - W,r) 
r V 9 

2(v g - v v ) v n 



+ iB j 



iB k 



(55) 



where erfc(z) is the complementary error function, the constant Cj k oc CjCl is fixed by the 
condition Ijj(L) = 1, and 

Ep-EAP) , , 

jK 1 . (56) 



B, 



Vj - 



The oscillation probabilities can now be obtained from eq. ( ]24j) . We will consider the limit 
of small a X7T , in which the result simplifies considerably. In the lowest non-trivial order in 
this parameter we obtain, after the proper normalization, 



Ijk(L) 



e 2P 



1 + 2.- 



v 9 Am jk 



l^hl V q - V n 2P 



For the z/„ survival probability we then find 



^ = e 4 + ^ + ^ 



[(cos 4> + £ sin (p) — Ar(£ cos (p — sin 0)] 



(57) 



(5? 
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2 Vn Am 2 h 

A _ 4 u _9 25 rr 

V27T - ^ 2P 



where 



The term independent of a X7T here coincides with the n p /v % 3> 1 limit of the z/ M survival 
probability (13TT) obtained in the case of neutrinos produced in decays of pointlike pions, 
whereas the term proportional to A v gives the correction due to the finite size of the pion 
wave packets. Note that eq. (1581) is actually the expansion of 

2cV 

Pw = c 4 + s 4 + [cos(0 - A n ) + £ sin(0 - Ar)] , (60) 

to the lowest non-trivial order in A n , which shows that to this order the effect of the finite- 
size pion wave packets reduces to an additional oscillation phase proportional to <r xn . 



4.2 Box- type shape factor of the pion wave packet 

Consider now a box- type shape factor for the pion wave packet: 

g n (x — v^t) = boxiv^t — x, d, 0) . (61) 

It corresponds to the pion wave packet of width d and the initial condition that the pion 
starts arriving in the decay tunnel at t = 0. From (H71) we have 

UEM -Ep) = i r • (62) 

Ej(p) -Ep + i^ 

To facilitate comparison with the Gaussian pion wave packet case, it will be convenient for 
us to express the parameter d through the effective spatial width of the pion wave packet 
a xn . The latter we always define as the coordinate dispersion in a given state, [(a; 2 ) — (x) 2 ] 1 ^ 2 . 
For the pion state flBTl) we thus find 

= £s ■ (63) 

Consider once again the limit Tl p /v n ^> 1. Substituting (1621) into fl52|) gives, in the 
leading non-trivial order in a X7T , 

e~ l 2P L r 1 v n Ami i 
1 - K L V 3 v g - v w 2P J 

The probability P MM is then again given by eq. ( |58|) but now with 

A -A - 1 v » ,t (RK\ 

V3 v g - v w 2P 
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Comparing eqs. ( 157)) and (164)) we see that they have the same structure: there is a term of 
the zeroth order in a X7T which coincides with the corresponding expression in the Yl p / tv ^> 1 
limit of the case of pointlike pions (see eq. ()28l) ). and the term which is linear in a X7T . The 
coefficients in front of a X7T in the cases of box-type and Gaussian type pion wave packets 
differ only by a numerical factor. The same applies to the expression for P MAt - in both 
cases it is given by the same eq. (1581) with the coefficients A n only differing by a numerical 
factor. Note that expressions ()57j) . (16%)) and (158)) were obtained in the limit of small A n and 
therefore cannot describe neutrino production decoherence caused by the finite size of the 
pion wave packet (which would correspond to A n > 1); they just describe the appearance of 
extra oscillatory terms in the transition and survival probabilities rather than suppression of 
such terms. We will consider the decoherence effects due to a xn ^ in section 16.21 For now, 
we will explain qualitatively the structure of the order a xn correction term in the oscillation 
probability (1581) and in particular its proportionality to v g /(v g — v n ). 



4.3 Small corrections due to a xn ^ 0: a qualitative analysis 

If the pion wave packet has a finite spatial width, there is a an additional contribution A<fi 
to the oscillation phase stemming from the fact that the neutrino production region (and 
the space-time interval of integration over the neutrino production coordinate) is different 
in this case. One has therefore to add A(f> to the oscillation phase <fi corresponding to the 
case of pointlike pions. 

Neutrino emission in a decay of an extended pion is schematically illustrated in fig. 2. 
The wave packet of the pion is represented by a band in the (t, x) plane, the muon wave 
packet is not shown. The slanted line represents the neutrino trajectory corresponding to 
the exact space and time localization of its detection. For neutrino emission from the space- 
time point A the additional oscillation phase as compared to neutrino emitted from the 
point B is given by the difference of the phases acquired by the mass eigenstates Uj and vj. 
along the segment AB shown in the figure: 

A0 = -[E 3 (P 3 ) - E k (P k ))At + {Pj - P k )Ax . (66) 

The quantities At and Ax are given by the projections of the segment AB on the t and x 
axes. Simple geometric considerations yield 

At = -^, Ax — (t X tt — ~ — ■ (67) 

Vg ~V V Vg~ V T 

Substituting (l6"7j) into ( 166)) and using eqs. (117)1 and ( ITS)) , we find 

v g - v n 2P 

This quantity does not depend on and therefore the same additional phase appears for 
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Figure 2: Schematical illustration of neutrino emission in pion decay in the case of finite- 
width pion wave packets. 

neutrinos emitted at any time. The quantity Ijk(L) can now be obtained from the one in the 
case of pointlike pions (eq. f l28|) ) by substituting </> + A<p. Taking the limit Tl p /v n ^> 1 
and keeping only the terms up to the first power in A0, for the probability P w we then 
obtain eq. ( 1581) with = —A(f> (cf. eqs. ( )58|) and (|60l)). The correction to the oscillation 
probability due to a xn ^ obtained here from simple qualitative considerations coincides, 
up to numerical factors of order one, with those found in the previous subsections for the 
Gaussian and box-type pion wave packets. The exact value of the numerical coefficient 
depends on the shape of the pion wave packet. 



5 Probability summation and finite-thickness targets 

We shall now study the case when a proton beam is incident on a target of finite thickness and 
will treat this case in the framework of incoherent probability summation. This procedure 
is just the one that is usually employed in the analyses of neutrino oscillation experiments 
in order to take into account the extended sizes of neutrino sources (and detectors). The 
results of this section will also turn out to be useful for comparing the probability summation 
approach with the coherent amplitude summation in the case of finite widths of the pion 
wave packets. We therefore examine this case in detail. 

Consider a bunch of protons of duration tt) Unc h with certain time distribution of protons 
within the bunch corresponding to a flux density F p (t). The protons are incident on a target 
of width D and density distribution n tar (x) which in general depends on the coordinate x 
along the direction of the beam. We assume that the protons produce pions on different 
scatterers incoherently. Consider the situation when the time of the pion production as well 
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as the time of neutrino detection are not measured and the neutrino signal is summed up 
over the bunch. The total number of pions produced in the layer (x, x + dx) of the target 
by the proton bunch is 

dM x) = K J dt A(x)F r (t) a n tar (x)dx = k „^> A<x M x )cl x , (69) 

where a is the pN— cross-section, A{x) is the factor describing the attenuation of the proton 
beam in the target by the time it reaches the point with the coordinate x, k n is the average 
number of pions produced in each interaction, and 

Fp nt = J dtF p (t) (70) 

is the integral proton flux. Notice that at this point the time dependence factorizes out 
and disappears from the whole picture: obviously, the obtained result does not depend on 
a specific time dependence of the proton flux. 

The number of muon neutrinos at the detector which originate from the decays of pions 
born in the layer (x, x + dx) of the target can be written as 

dn n (x) 1" ' dx'e-^ {x '- x) Pf™ A {L-x'). (71) 

J X 

Here x' is the distance from the front end of the target (from the side of the incoming proton 
beam) to the point of the pion decay, L is the distance from the target's front end to the 
detector, and P^ and (L — x) is given by eqs. ()42l) and (132]) with L replaced by L — x. The 
total number of muon neutrinos at the detector site, which originate from the decays of 
pions produced over the whole target, is 

N v = aF; nt K £ dz n tar .{z)A{z) j'* dxe~^ x ~ z) P^ nd (L - x) . (72) 
Let us consider some special limits of this expression. 



1. Very thin target: 

n tar (z)=n 5(z), A(0) = 1. (73) 

In this case ^ 

N„ = F; n %n [ " dxe-^ x P^ nd (L - x) , (74) 
Jo 

and the integration leads to the effective oscillation probability that coincides with the 
one in eq. (JSJ (see (|35|)). 

2. Short decay tunnel: l p v n /T. In this case 

N v w F; n % [ D dz n tar {z)A{z) f* dxP s ^ A {L - x) . (75) 

J0 J z 
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Assuming for simplicity that n tar (z) = n tar = const and A(z) = A Q = const, we find 



N v « F; nt Kn tar A I dz I dxP s ;™ d (L - x) . (76) 



o 



Straightforward integration yields 
N, 



D , . , 7iD 



sin — sin(0 — <p p ) + — (sin — sin(0 — <p p )) — - — cos < 



v — T — vr 



(77) 

where we have retained only terms up to the first order in D. In the limit D — > this 
expression reduces to the one in eq. (jUJ). 

3. Extended target and well localized protons. Assume that the decay tunnel is long, 
so that /decay = v n /T l p . This is the same limit as the one used in the coherent amplitude 
summation calculations with finite- width pion wave packets carried out in sections 14.11 and 
14.21 We also assume that the density of the target is constant and that there is no absorption, 
i.e. n tar (z)A(z) = const. Thus, we have the setup which corresponds to the one that in the 
amplitude summation approach led to the box-type pion wave packets. From eq. ( 172|) we 
find 

N„<x I dz I dxe-^ {x ~ z) P^ nd (L - x) . (78) 

The integration yields 



o 



N„ A . a . 2c 2 s 



2„2 r A™2 



PeS = « C + S + 



/\rn 

cos0 + ^sin0 ^-D(^cos0 — sin0) , (79) 



where only corrections of the first order in D are kept. 

Now let us compare the above results with those of the coherent amplitude summation 
approach with finite spatial width of the pion wave packets. 

What could be an analogue of the finite-width pion wave packet in the case when one 
performs the summation of contributions of neutrino production at different points at the 
level of probabilities? In that case the standard oscillation probability is used, and the notion 
of the pion wave packet does not apply. Nonetheless, one can introduce an analogy between 
the incoherent probability summation approach and the coherent amplitude summation 
with finite a xn in the following way. 

In sec. HI when considering the oscillations of neutrinos produced in pion decays in the 
amplitude summation approach, we estimated the width of the pion wave packets assuming 
that the pions are produced by proton collisions with nuclei of a target and that these 
nuclei are well localized and are essentially pointlike, whereas the protons are described by 
extended wave packets (see eqs. (Q and (j7])). However, one can instead imagine that the 
protons are pointlike, whereas the nuclei of the target are localized in relatively large spatial 
region of size D, so that their wave packets are characterized by the spatial width a X N ~ D. 
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In this case eqs. (j3J), (jlj) and (jSJ) yield <r x p = cr xp ~ 0, fp ~ v p , <r evr ~ v p /D, and from eq. 
we find, instead of (|7j), 



D . (80) 



On the other hand, in the probability summation approach one can identify the spatial size 
of the localization region of the target nuclei with the extension of the target in the direction 
of the proton beam. Therefore, when one employs the probability summation approach, the 
role of an analogue of the finite-size pion wave packets could be played by the finite-thickness 
target, with the correspondence between the thickness of the target D in the probability 
summation approach and the width of the pion wave packet a xn in the amplitude summation 
one established by eq. (I80p . Such a calculation was carried out in example 3 above, and the 
result is given in eq. (1791) . Note that the obtained expression for the effective muon neutrino 
survival probability P e ff has the same structure as the probability P^ in eq. (1581) . which 
was derived in the coherent amplitude summation approach in the leading non-trivial order 
in a X7T . The two expressions, however, differ by the coefficients in front of the second terms 
in the square brackets: while in the amplitude summation framework with box-type pion 
wave packets it was given by A whox of eq. (16"51) . in the probability summation approach of 
this subsection it is 

Am 2 Am 2 v p 
4P AP v p — v n 

Here in the last equality we have used relation ( 180]) which establishes the correspondence 
between the proton target thickness D of the probability summation approach and the 
width of the pion wave packet of the amplitude summation one. It is interesting to note 
that the expressions in eqs. ( l65l) and (IHTl) differ only by the replacement of the neutrino 
group velocity v g by the proton velocity v p and by the numerical factor a/3/2. 

The fact that eqs. (1551) and (179")) have similar structure is actually not surprising: it can 
be shown under very general assumptions that the corrections due to the finite thickness 
of the proton target should have such a structure, irrespective of whether the coherent 
amplitude summation procedure or the incoherent probability summation one is employed. 
It should be stressed, however, that the coefficients of the correction terms are different. As 
will be shown in section 16.11 this is a reflection of the fact that the equivalence between the 
coherent and incoherent summation approaches in general only holds in the case of pointlike 
pions and localized in space and time neutrino detection. It can also be shown that the 
equivalence holds even in the case of extended pion wave packets provided that both the 
neutrino and muon detection processes are perfectly localized (see section [9]). 
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6 An alternative approach to calculation of oscillation 
probabilities. Finite-width pion wave packets 



We shall now develop a different calculational approach to the coherent amplitude summa- 
tion procedure, which is based on a different order of integrations involved in the calculation 
of Ijk(L) and is valid for arbitrary shapes of the pion wave packets. This approach will be 
useful for explaining the equivalence of the results of the amplitude and probability summa- 
tions of the contributions of neutrino production at different points, which was established 
above for the case of pointlike pions and perfectly localized neutrino detection. It will also 
be helpful for studying the cases when the muon emitted alongside the neutrino in pion 
decay is detected and when the interactions of pions in the bunch between themselves or 
with other particles inside the neutrino source are taken into account. 



6.1 General formalism 

We start with eq. fl50l) for Ijk{L) and substitute into it expression ( [12]) for the momentum- 
space neutrino wave packet. The pion wave function is taken in the form (j45p with an 
arbitrary shape factor g n (x — v n t). The resulting expression contains integrations over two 
coordinates, two times and a momentum. Integrating first over the momentum, we obtain 

Am 2 , r r 

I jk (L) = C e-' l ^ L / dxtdh / dx 2 dt 2 g n (x x - v n h) g*(x 2 - v j 2 yl^(P)-Ep]ti~i[E k (P)-E P }t 2 

Am 2 , „ 

x e i^[(x,+X2)-v g {t 1+ t 2 )]~^t 1+ t 2 ) ^ _ x ^ _ ^ _ ( g2 ) 

Here Co is a constant, and the integrations over x\ and x 2 are performed in the interval 
[0, lp\. The integrations over time are carried out over the whole interval (— oo, oo), which is 
justified if g n (z) decreases sufficiently rapidly for large \z\ (see the discussion after eq. (H7|) ). 
After the integration over the time variables, eq. (1821) becomes 



. Am 2 fc t r r 
I jk {L) = Ce 1 -*p J dxi j dx 2 e 



Am 2 , „ 

xj^-^+^-^^+^G^xt - x 2 ) , (83) 

where C is a constant and the effective shape factor G n (xi — x 2 ) is a convolution of two g^ 
functions (taken at shifted arguments) with an exponential factor: 

/qj qj qj qj Yy 

dygn(y + 9 2v w (xi - x 2 )) g*(y - g ^ n (xi - x 2 )) . (84) 

The function G^{x\ — x 2 ) satisfies 

G;(z) = G n (-z) . (85) 
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If is real, then G n is also real and is an even function of its argument. 

Note that the shift of the arguments of the two g m functions in (|84|) is proportional 
to (v g — v n )/v g . The presence of the two g n functions with shifted arguments in (154"]) is 
a reflection of the fact that we have a coherent summation of the amplitudes of neutrino 
production at different points here (there are interference terms between g n and g* which for 
the same value of y have different arguments). For pointlike pions g n (x — v^t) oc S(x — v n t), 
and instead of the product g n (zi)g*(z2) we get \g n ( z i) \ 2 , meaning that there is no interference 
terms and therefore the amplitude summation reduces to the probability summation. This 
gives the explanation of the equivalence of the two summation approaches, established in 
sections |2] and [3j We will discuss this point in more detail in the Discussion section. 

If the pion shape factor g n (x — v n t) is characterized by a width a X7T (i.e. g n (x — v^t) 
quickly decreases when \x — v n t\ becomes large compared to cr X7r ), then, as follows from its 
definition, the function G n (xi — 22) is typically centered around the zero of its argument 
and is characterized by the width [v g /(v g — v n )}a xv . 

Equation can be transformed into the form0 

rip pip Am2 fc 

I jk {L)=C dzK jk (z) die"^ M -- (l - 2) , (86) 

JO Jz 

where 

K jk (z) = [G„(z)e ilEk(P) - Ep] % + Gl{z)e- l[E > {P) - Ep] ^]e-^ . (87) 

Eqs. (186]) and (IHTj) are the final result of the calculational approach presented in this section. 
As usual, the constant C in ( 186]) should be determined from the condition Ijj(L) = 1. Note 
that, even though formally the integration over z in (186]) extends from to l p , in reality it 
is effectively limited from above by z max ~ [v g /(v g — v n )]cr X7T due to the properties of G n (z) 
(provided, of course, that z max < l p ). For the special cases of Gaussian and box-type pion 
wave packets and Tl p /v n 3> 1 eqs. (181]) - (I8T]) reproduce the results obtained in sections H~T1 
and O 

It is interesting to note that there is a striking similarity between the general result (1861) 
for Ijk{L) in the coherent amplitude summation approach and eq. (I72j) obtained for finite- 
thickness targets in the incoherent probability summation one. There would be a complete 
correspondence between the two results if the quantity Kj k (z) in (I8B]) did not depend on 
the indices j and k: Kj k (z) — > K(z). In that case it would be possible to multiply the inner 
integral in (186]) by U*jU a kU/3jU* k and perform the summation over j and k, which would 
produce P^ nd (L — x) in the integrand. The same operation applied to the left-hand side 
of (186]) yields, according to (124]) . the oscillation probability P a/ j(L). The correspondence 
between the two results would then be exact provided that one identifies n tar (z)A(z) of 

7 For this, one has to go to the integration variables z = x\ — xi and x = Xl + X2 ; write the outer integral 
over z as dz — dz + J Q ' P dz, change the integration variable in the first of these integrals according 
to z — > — z, and shift the integration variable in the inner integral over x according to x — > x + z/2. 
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eq. ( !72|) with K(z) of eq. (1861) . The quantity Kjt(z) is independent of the indices j and 
fc, e.g., in the case of pointlike pions. Indeed, from ( |84|) and (!H7|) we find that in this case 
Gt T (z) oc 5(z) and Kjk(z) = K{z) = 2G v {z). As we already know, for pointlike pions the 
results of the amplitude and probability summation approaches indeed coincide. 



6.2 Large a X7T and production decoherence 

If the spatial width of the pion wave packets a xn is sufficiently large, one can expect produc- 
tion decoherence due to the averaging of the oscillation phase along the pion wave packet. 
In analogy with the parameters £ and <j) p in eq. fl39|) . the decoherence parameter in this case 
is expected to be 

Am 2 . Ax , . 

&r = -^p- Ax = 271— , (88) 

^osc 

where the quantity Ax defined in eq. f)67p is the spatial size of the region over which the 
additional averaging of the oscillation phase occurs. Note that the parameter <p p coincides 
with the introduced earlier quantity A n up to a numerical factor of order one. Neutrino 
production decoherence due to large a xn might be expected to manifest itself through the 
suppression of the oscillatory terms of the flavour transition probabilities by ^-dependent 
factors that vanish in the limit <f> % — > oo. However, large-cr X7r expansions of the oscil- 
lation probabilities in the cases of Gaussian and box-type pion wave packets, for which 
we have closed-form expressions, do not reveal any suppression of the oscillations due to 
(^■-dependent factors. How can this be understood? 

Let us first consider the limit Tl p /v w ^> 1. As follows from eq. (186|) . Ijk(L) oc 1/(1 
in this case. The oscillatory terms of the probabilities P a p{L) will then be in any case 
suppressed if £ 3> 1. Let us now notice that £ <f> % , which follows from (fTTTm It should 
be stressed that (TTTj) is a necessary condition for not including the pion production process 
in the description of neutrino oscillations, which is assumed to be satisfied throughout this 
paper. Thus, in the limit of large 4> p we have £ ^> (f) p ^> 1, and the oscillations are suppressed. 

Similarly, in the opposite limit Tl p /v n 1 the oscillations are suppressed provided 
that (p p ^> 1 (see the discussion in section 12. 5j) . Note that (f> p ^> (f> n because l p ^> a X7T . 
Therefore, in the limit Ylp/v^ 1 the oscillations will be suppressed for large 0^ because 
of (pp ^> <p n 1. It is easy to see that for large (f> n the oscillations are also suppressed in 
the case Tl p /v n ~ 1. We conclude that in the case <p n ^> 1 the oscillations are quenched as 
expected, even though not directly by 0„--dependent factors. 



8 Except for pions of extremely high energies, which are currently inaccessible. Indeed, £,/4>n = 
(Jdec&y / &x-k) [(v g ~v^)/v g ]. The first factor on the right hand side of this equality is extremely large because 
the pion decay length /decay = v n /F is a macroscopic quantity, whereas the width of the pion wave packet 
axn is microscopic. The condition £/4>tt < 1 therefore can only be satisfied for extremely small values of 
V g — v-k « 1 — Wtt- From a similar argument it also follows that (f> p 3> (j> v . 
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7 Effects of muon detection and pion collisions 



Up to now we have been considering only the situations when the muon produced together 
with the neutrino in the pion decay is neither directly detected nor interacts with the 
medium, and the pions in the bunch do not interact with each other or with other particles 
which may be present in the neutrino source. We shall now lift these restrictions. 

As in section |2} we shall be assuming here that the parent pions can be considered as 
pointlike and are described by the wave function ([9]). Consider the case when the muon 
produced alongside the neutrino in the pion decay is "measured", either directly by a dedi- 
cated detector or through its interaction with the particles of the medium. In this case the 
muon is localized by the interaction and therefore it should be described by a wave packet 
rather than by a plane wave. We adopt the wave function of muons of the form (jSJ): 

^(x, t) = e iKx - iE ^ g,[(x - x s ) - v,(t - t s )} . (89) 

Here g^[{x — xs)—v^{t—ts)\ is the shape factor of the muon wave packet, which is determined 
by the muon detection process. It is assumed to have a peak at the zero of its argument 
and decrease rapidly when the modulus of the argument becomes large compared to a XfJi , 
where a xil is the spatial width of the muon wave packet. Our choice of the argument of 
g M corresponds to the initial condition that at the time t = t$ the peak of the muon wave 
packet is at x = x$- We shall choose x$ to be the coordinate of the center of the muon wave 
packet at the neutrino production time. Since we assume the pions to be pointlike and the 
point xs should lie on the pion trajectory, xs and ts must be related through xs = v^ts- 
Obviously, xs <l p - 

Repeating essentially the same calculations as in section El we arrive at a very simple 
and compact expression for Ijk(L): 

I jk (L) = C dx\g,((x-xs)-^^)\ 2 e- l ^ {L - x) - r ^. (90) 

The argument of the shape factor function g^ here implies that the effective width of the 
muon wave packet that enters into eq. (190]) is actually 



(Txfj, = <?xn ~ • (91) 



Note that this expression has a structure analogous to that of Ax in eq. fl67j) and allows a 
simple geometric interpretation in terms of the plots similar to the one in fig. [2j 

The calculation leading to (190]) was carried out in the quantum-mechanical approach 
with coherent amplitude summation of the contributions of neutrino production at different 
points along the pion path. Nevertheless, in the course of the calculation one had to integrate 
an expression containing \g n (z) | 2 rather than g^(zi)g*(z2). As was stressed in section |6j 
this reflects the fact that for pointlike pions and localized neutrino detection the results 
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of the amplitude and probability summation approaches coincide. Thus, the equivalence 
of these two approaches in this case holds true irrespective of whether the charged lepton 
accompanying the neutrino production is detected or not. This is also confirmed by the 
presence of the squared modulus of in the integrand of eq. ( 1901) . We will discuss how the 
charged lepton detection can affect the equivalence of the two approaches in the cases when 
either the pion is not pointlike or the neutrino detection process is delocalized in section [91 

Let us now study the result in eq. fl90|) in several limiting cases. First, we note that 
in the limit a Xfl — » oo, which corresponds to the plane- wave approximation for the muon 
employed in sections [2] - [61 the shape factor of the muon wave packet — > const, and we 
recover the previously found expression ([28]) for Ijk(L) and therefore the previous results 
for the oscillation probabilities. 

Consider now the case a XIM -/)■ oo. Eq. (|9"U|) then actually describes oscillations of a 
"tagged" neutrino, i.e. of a neutrino produced together with the muon which was detected 
and whose production coordinate was found to be xs with the accuracy given by a Xfl . Let 
us first consider the limit a x/1 — > 0. In this case in the integrand of eq. (I90I) the function 
g M oc 5(x — xs), i.e. the muon detection exactly localizes the neutrino production point. 
Eq. (J9QD then yields 

I jk (L) = const. e" r ^ e -i^(L-x s ) _ ( Q2 ) 

The real exponential factor here just describes the depletion of the pion flux by the time they 
reach the point x = x,s', this factor can be absorbed in the overall normalization of Ijk{L) 
which does not affect the transition probability. The usual normalization procedure yields 

Ijk{L) = exp[— % ™p k (L — xs)], which leads to the standard probability P^ and (L — x s ) that 
describes v a — > z/g oscillations between the neutrino production point xs and the detection 
point L. The production decoherence effects are absent in this case. 

Consider now the case when the effective spatial width of the muon wave packets a xlx 
is neither infinite nor vanishingly small. For definiteness, we shall consider Gaussian muon 
wave packets with the shape factor g^(x — v^t) = exp[— ^~J 2 M ^ ], which gives 

9Mx - xs)^—^) = e . (93) 

Let us now consider Ijk(L), as given by eq. ( 1901) with g^ from (1931) . in some limiting 
cases of interest. First, for very large a Xfl one expects to recover the results of section [21 
which were obtained for the plane-wave description of the muon, i.e. in the limit a x/1 — > oo. 
Indeed, taking into account that < xs < l p , its is easy to see that for a XfM 3> l p the function 
g M in the integrand of ( 19"U|) can to a very good accuracy be replaced by unity, so that we get 

rip Am ?fc 

I jk (L) = C dxe- 1 -^'-^, (94) 
J o 

which is just the result found in section [5] where the muon was assumed to go "unmeasured" 
and described by a plane wave. Indeed, direct integration shows that f )9"4"|) leads to (T2"Bj) . 
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Eq. (194)) is also obtained in the cases when a x/J/ < l p , provided that the main contribution 
to the integral in (190]) comes from a relatively small region x < x c with x c <C l p , and in 
addition a XfJi 3> x c . For instance, if Tl p /v n ^> 1, the integrand of (J90l) is strongly suppressed 
for x ^> v n /T due to the exponential decay factor. Thus, for T/v n > Am 2 k /2P we have 
x c ~ ^/r. In this case the factor ^ in the integrand can be replaced by unity and the 
results of section [2] are recovered provided that 

( v (v \ ~i 
cr XjU > max|— , ^— x 5 J J. (95) 

Analogously, if (Am 2 jk /2P)l p > 1 and Am 2 k /2P > T/v n we have x c ~ 2P/Am 2 k (for 
larger values of x the integrand of ( 1901) is fast oscillating and the corresponding contributions 
to the integral are strongly suppressed). In this case the results of section [2] are recovered 
provided that 

( 2P (IP \ 1/21 . . 

}• (96) 



cr^ > max 



Let us now consider the case of small a XfJL . In the limit 



a Xfl < min{/ p , x s } , (97) 

one can extend the integration in (190|) over the whole infinite axis of x. The integration 
with from ( 19~T|) then gives, after the usual normalization procedure, 



Am ?fe ,, r 9 n l f Am ?fc\ 2 -9 



Let us first look at the exponential phase factor in this formula. In the absence of the 
second (real) exponential factor, it would just describe the standard v a — >■ up oscillations 
between the points with the coordinates xs — (T/v^a 2 and L. The term (T /v^a 2 ^ in the 
coordinate of the initial point describes the shift of the position of the center of the neutrino 
production region due to the pion instability. As follows from ( I98p . its contribution to the 
oscillation phase can be neglected provided that 

\ V 2 x 1/2 

2P vA \ x - 



(7 ^ < l Am%T J ~ [ 2 n loscldccay ) ' (99) 

The second exponential factor (1981) leads to the suppression of the oscillations in the case 
(Am 2 k /2P)cr XIM ^> 1, i.e. it describes possible production decoherence effects. In particular, 
for the muon neutrino survival probability in the 2-flavour case we obtain, assuming that 
(J99D holds, 

P m = c 4 + s 4 + 2s 2 c 2 e -K^) a «, cos {^f(L - Xg j) . (100) 
Thus, in the limit (19~T1) the decoherence parameter is ^p-o^. 
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Finally, consider the case when the muon interacts with the medium but there are no 
muon detectors and the muon position is not measured. Equivalently, one can imagine that 
the muon position is measured but the correspondence between the detected neutrino and 
the coordinate of the production point of the accompanying muon is not established. In 
both these cases the neutrinos are not tagged, and one has to integrate ( )90l) over xs- The 
integration of the squared modulus of in f )90|) over x s just gives the normalization constant 
of this function which does not influence the oscillation probabilities, and we simply recover 
the results obtained in the case when the muon is not detected. 

From the above discussion of the muon detection case it is clear what happens if the 
interaction of the pions in the bunch between themselves or with other particles which may 
be present in the neutrino source is taken into account. Assume first that this interaction 
identifies the individual pion whose decay produces a given neutrino. For example, the pion 
decay may lead to some recoil of the neighbouring particles which may be detected. This 
would localize the coordinate of the neutrino production point within an uncertainty of order 
of the inter-pionic distance (or, correspondingly, the distance between the pion and the other 
neighbouring particles in the source) r , and would lead to neutrino tagging. The production 
decoherence parameter in this case is (Am 2 /2P)r . If, however, the information about the 
interaction between the decaying pion and the surrounding particles is not recorded and not 
used for neutrino tagging, the oscillations occur in exactly the same way as if pions did not 
interact with each other or with other particles. 

8 Implications for neutrino oscillation experiments 

We shall now estimate possible effects of production decoherence for a number of past, 
ongoing and forthcoming or proposed neutrino oscillation experiments. 

As was pointed out in section I2.1[ in realistic situations the spatial widths of the wave 
packets of parent pions are much smaller than all the length parameters of physical interest in 
the problem, therefore we will be using here the formulas from section[2]obtained for pointlike 
pions. It was demonstrated in that section that in the case Yl v jv v 3> 1, when most of the 
pions decay before reaching the end of the decay tunnel, the production decoherence effects 
are governed by the parameter £ defined in eq. ( 133]) (see also [3]), whereas in the opposite 
case Tl p /v n <C 1 the decoherence parameter is P defined in (1321) . In table 1 we give the 
parameters Tl p /v n , <p p and £, along with some other relevant physical characteristics, for a 
number of experiments [SllIllT5llISlinilIHllI2ll2nil2Ill22ll^. If not otherwise specified, 
we assume Am 2 = 2 eV 2 for all experiments. For a number of experiments (NOMAD, 
CCFR, CDHS) we adopt the values of Am 2 that correspond to the maximal sensitivity of 
these experiments. For /3-beams we consider a short-baseline setup with L = 130 km, the 
neutrino energy in the ion rest frame Eq = 2 MeV, ion lifetime tq = 1 s and 7 = 100. 

A few comments are in order. Obviously, the results we have obtained for the case when 
neutrinos are produced in pion decays are also valid for decays of any other parent particles, 
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Table 1: Production coherence for a number of experiments [12-22]. Unless otherwise spec- 
ified, Am 2 = 2 eV 2 . For /3-beams we adopt E = 2 MeV, r = 1 s, 7 = 100. 



Experiment 


(£„>(MeV) 


L(m) 


l p (m) 


Zdec(m) 


^osc(m) 


Tlp/vp 


(t> P 




LSND 


~40 


30 








50 


- 








KARMEN 


~40 


17.7 








50 


- 








MiniBooNE 


~800 


541 


50 


89 


992 


0.56 


0.32 


0.56 


NOMAD 


2.7- 10 4 


770 


290 


3009 


33480 


0.1 


0.054 


0.56 


(20 eV 2 ) 
CCFR(10 2 eV 2 ) 










3348 


0.1 


0.54 


5.64 


5-10 4 


891 


352 


5570 


1240 


0.06 


1.78 


28.2 


CDHS 


3000 


130 


52 


334 


3720 


0.155 


0.088 


0.56 


(20 eV 2 ) 










372 


0.155 


0.878 


5.64 


K2K 


1500 


300 


200 


167 


1861 


1.2 


0.68 


0.56 


T2K 


600 


280 


96 


66.4 


744 


1.45 


0.81 


0.56 


MINOS 


3300 


1040 


675 


368 


4092 


1.84 


1.04 


0.56 


NOz/A 


2000 


1040 


675 


223 


2480 


3.03 


1.71 


0.56 


/3-beams 


400 


1.3-10 5 


2500 


3-10 10 


496 


8.3-10" 8 


31.7 


3.8-10 8 



including muon decay, provided that the particles produced alongside the neutrino are not 
detected. Next, we note that for LSND and KARMEN experiments, in which neutrinos 
were produced in decays of parent muons at rest, the sizes of the neutrino sources are very 
small compared to the baselines and oscillation lengths and therefore can be considered to 
be essentially zero. This means that the produced neutrino states are characterized by large 
energy and momentum uncertainties, i.e. for these experiments the production coherence 
condition ([[]) is therefore satisfied very well. Notice that the parameter £ is practically 
energy independent (it depends essentially only on the neutrino production process and 
Am 2 ). This is because the decay width in the laboratory frame V oc 1/E, and in the 
denominator of (133]) its energy dependence cancels with that of P ~ E. The only remaining 
energy dependence is through the pion velocity v n (or in general through the velocity of the 
parent particle vp), which is very weak for relativistic parent particles. For it — > \w decay in 
flight and Am 2 ~ 2 eV 2 the parameter £ is always ~ 0{1). This does not necessarily mean 
that production coherence is always violated in this case. Indeed, as discussed in section 
12. 5[ for l p < /decay production decoherence is governed by <fi p rather than by £. 

As can be seen from the table, production decoherence effects should be noticeable for 
MiniBooNE, NOMAD (Am 2 = 20 eV 2 ), CCFR (Am 2 = 100 eV 2 ), CDHS (Am 2 = 20 
eV 2 ), K2K, T2K, MINOS and NOz/A. Very large decoherence effects are expected for short- 
baseline /3-beams with the parameters quoted in the table. For illustration, in fig. [3]we show 
the oscillation probabilities 1 — P MAt for the T2K and CCFR experiments for the parameters 
indicated in the table and = 0.2. The red curves show the actual oscillation probabilities, 
whereas the blue curves were obtained for neutrino emission from a single point located in 
the middle of the decay tunnel, i.e. correspond to the probabilities obtained neglecting the 
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Figure 3: Probabilities 1— P MM for the T2K experiment (left panel) and the CCFR experiment 
(right panel) for parameters indicated in table 1 and = 0.2. Red curves: the actual 
oscillation probabilities, blue curves: the oscillation probabilities in the case of neutrino 
emission from a single point located in the middle of the decay tunnel. Neutrino energies 
are in MeV. 

decoherence effects. Another example of production decoherence can be found in fig. 3 of 
ref. |24j, where v e — > v s oscillations in low-energy /3-beam experiment with 7 = 30, l p = 10 
m and L = 50 m were considered. Effects of suppression of the oscillations due to the 
averaging of the neutrino production coordinate over the straight section of the storage ring 
and of the neutrino absorption coordinate along the detector can be clearly seen. 

9 Discussion 

In the present paper we have studied in detail neutrino production coherence, which ensures 
that the emitted neutrino state is a coherent superposition of different mass eigenstates, as 
well as its implications for neutrino oscillations. We considered neutrino production using 
71 — > fiu decay as an example. However, our analysis also applies to neutrinos produced in 
any other decay process (such as A-decay, /i-decay, /3-decay, etc.), provided that particles 
accompanying neutrino production are undetected, as it is usually the case. 

We studied neutrino production coherence in two completely different approaches. In 
one of them, based on the quantum-mechanical wave packet formalism, we found the tran- 
sition amplitude A a /3(L,t) by summing the amplitudes of neutrino production at different 
points along the path of the parent particle. The oscillation probability P a p(L) is then 
found by integrating the squared modulus of the obtained amplitude over time. In the 
second approach we assumed that each individual neutrino production event is fully coher- 
ent, so that the evolution of the produced neutrino should be described by the standard 
oscillation probability. However, the neutrino emission can occur (with some probability) 
at any point within the neutrino source, and therefore the effective oscillation probability is 
obtained by the proper averaging of the standard probability over the neutrino source. In 
both approaches we assumed that the neutrino detection process is perfectly localized. 
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We have found that in the case when the parent particles, in decays of which the neutri- 
nos are produced, can be considered as pointlike, the results of the two approaches exactly 
coincide. This happens despite the fact that the first approach is fully quantum mechanical, 
while the second one, based on the probability summation, is essentially classical. We thus 
confirm the conclusion of ref. [3] where, along with the probability summation, a simplified 
treatment of the quantum mechanical coherent amplitude summation was given. 

It should be stressed that the equivalence of the two approaches studied in this paper is 
different from the one discussed in [5], [25]. In those papers the authors have pointed out that 
it is impossible to distinguish experimentally between an ensemble of neutrinos described by 
wave packets, each with the energy distribution amplitude g{E), and a beam of plane- wave 
neutrinos, each with well defined energy, but with the energy spectrum of the beam given 
by the squared modulus of the same function g(E). In our case we establish the equivalence 
between a beam of neutrinos described by wave packets and a beam of pointlike neutrinos, 
which are just the opposites of the plane waves. 

How can one understand the equivalence of the two approaches studied in this paper 
which takes place in the case of pointlike pions? Although the position of the pion decay 
(and neutrino production) point is not a priori exactly known, having uncertainty of order 
l p for Tl p /v n <C 1 and /decay = v n /T for Tl p /v n ^> 1, the spatial size of the production region 
for each individual production event is very small: it is given by the size of the smallest 
wave packet of the particles participating in neutrino production, in our case of the pion. 
For pointlike pions the space-time localization of the detection process actually allows one, 
for each detection event, to pinpoint the coordinate of the neutrino emission. Thus, in this 
case there is no quantum mechanical uncertainty in the coordinate of neutrino production 
and therefore no interference between the amplitudes of neutrino emission from different 
(even closely located) points. This explains why our results obtained through the coherent 
summation of the amplitudes of neutrino production along the neutrino source coincide with 
those found by the simple incoherent summation of the oscillation probabilities. 

From the above argument it follows that one can expect some deviations between the 
results of the coherent amplitude summation and incoherent probability summation ap- 
proaches if pions are described by wave packets of finite size o~ X7T . This is indeed confirmed 
by our treatment of the finite o~ X7T case in sections H] and [6j We have found that for small 
Ar ~ [ v g/( v g ~ v ir)](Am 2 /2P)a xn the oscillation probabilities get extra oscillatory terms 
proportional to A n . For different shapes of the pion wave packets these corrections differ 
only by numerical factors (cf. eqs. fl59l) and fl65l) ). At the same time, for ^> 1 the 
neutrino production coherence is violated, and the oscillatory terms in the oscillation prob- 
abilities are strongly suppressed. From the above discussion it should be rather obvious that 
the amplitude summation and probability summation approaches should lead to different 
results also in the case when the detection process is not perfectly localized, i.e. when the 
particles participating in neutrino detection are described by wave packets of finite size. 
This should be the case even if the production process is perfectly localized. 
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The above points as well as the role of detection of the charged lepton produced alongside 
the neutrino are illustrated by fig. HJ In this figure we present the space-time diagrams that 
correspond to six different experimental setups in the case of neutrinos produced in it — > \w 
decays. The size of the region of coherent amplitude summation in each case is determined 
by the interplay of three factors, namely, whether the pion can be considered as pointlike or 
not and whether the neutrino and muon detection regions can be considered to be space-time 
points or notE Whenever this region degenerates to a point, the coherent and incoherent 
summation approaches yield identical results. For simplicity, in the cases when the neutrino 
and the muon detection processes are not fully localized, we display them as being extended 
in time but localized in space. More general situations where these processes are delocalized 
both in space and time can be readily studied; however, the corresponding results will not 
modify our qualitative conclusions. The six panels in fig. H] thus illustrate the following 
situations: 

a) The pion and the neutrino detection region are both pointlike, the muon goes unde- 
tected. The first two conditions are sufficient to identify the point where the neutrino 
was produced, thus eliminating any quantum mechanical uncertainty in the emission 
coordinate. The oscillation probability found through the coherent amplitude sum- 
mation must be identical to the one obtained by incoherent probability summation. 

b) The pion and the neutrino detection region are both pointlike, the muon detection 
region has a finite size. The latter does not restore coherence of neutrino emission 
from different points (illustrated by the dotted lines in the figure). The corresponding 
contributions have to be summed at the level of probabilities. 

c) The pion is of finite extension and so is the muon detection region, whereas the neu- 
trino detection is pointlike. A one-dimensional region AB is formed as the intersection 
of the regions corresponding to these three conditions. It is impossible in principle 
to determine from which point in the segment AB the neutrino was emitted, and 
the amplitudes of neutrino emission from all such points thus interfere. The ampli- 
tude summation produces the results that are different from those found through the 
probability summation. 

Depending on the degree of its derealization, the muon detection process may reduce 
the length of the segment AB compared to the case when the muon is undetected 
(cf. figs. [2] and Hb) and thus diminish the effects of averaging of neutrino oscillations 
caused by neutrino production decoherence. 

d) The pion is of finite extension, the neutrino detection region is also finite, but the 
muon detection region is pointlike. As in c), the region of amplitude summation is a 
segment, and there is no equivalence with incoherent probability summation. 

9 Clearly, this can be made more precise. For the pion, one needs to compare Ox^Vgjivg — v v ) with / osc . 
For the detection processes, it is necessary to compare the size of the detection region with ^ osc . 
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e) Both pion and neutrino detection regions are of finite size, while the muon goes un- 
detected. A 2-dimensional region of amplitude summation is formed; the amplitudes 
of neutrino emission from different points of this region interfere. 

f) Same as in e), but with muon detection (finite detection region). As in case e), the 
region of amplitude summation is 2-dimensional. Its size may be smaller than in case 
e) due to constraints from the muon detection. 

Not shown in fig. H] is one more case when the results of the amplitude summation and 
probability summation approaches coincide: neutrino detection region is of finite extension, 
whereas the pion is pointlike and the muon detection is perfectly localized in space an time. 

Are deviations between the results of the coherent amplitude summation and incoherent 
probability summation approaches experimentally observable? Consider the case of perfectly 
localized neutrino detection, non-zero spatial widths of the pion wave packets and no muon 
detection (fig. [2]). The corrections to the oscillation probabilities due to a xn ^ are governed 
in this case by the parameter A n which, as follows from ( 1591) or ( |65i) . can become sizeable 
for extremely high energies of the parent pion. This requires 



For instance, if the width of the proton wave packets a xp ~ 10~ 4 cm, one can expect also 
i$ 10~ 4 cm (see eq. (I7j)); for Am 2 ~ 1 eV 2 the parameter A n will then become of order 
one for pion energies E„ > 10 3 TeV. Such energies are not feasible, and even if they were, 
the corresponding neutrino oscillation lengths would be far too large for any oscillation 
experiment. Another possibility to make A n sizeable would be to increase significantly the 
spatial width of the wave packets of ancestor protons, which would in turn increase the 
values of a X7r . However, it is not clear how this can be achieved Jjjjl 

The two approaches to production decoherence effects that we followed in this paper 
describe the same physical phenomenon - the suppression of the oscillating terms in neutrino 
transition and survival probabilities due to derealization of the production process. The 
nature of this suppression, however, is different in different situations. 

In general, the proper description of the decoherence phenomenon requires using the 
quantum-mechanical approach with coherent summation of the amplitudes of neutrino pro- 
duction at different points. In this framework one considers the violation of production 
coherence due to lack of localization of each individual neutrino production event. This 
derealization of neutrino production is related to the fact that the exact coordinate of the 
neutrino emission point is not known. The decoherence parameter is therefore essentially 
the ratio of the size of the neutrino production region and the neutrino oscillation length. As 
was shown in section [21 in the limit Tlp/v^ <C 1 the decoherence parameter is <p p = 2ttI p /1 osc , 

10 It might still be possible, however, to achieve relatively large widths of the wave packets of the parent 
particles if neutrinos are produced in muon decays. 
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Figure 4: Space-time diagrams corresponding to six experimental setups. The regions of 
coherent production in each case are shown in deep red (see text). 
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whereas for Tl p /v n 3> 1 it is £ = 27r/deca y //osc- The decoherence parameters can also be 
represented as the ratio of the energy difference of different neutrino mass eigenstates and 
the energy uncertainty inherent in the neutrino production (see section |2~5T) : AE/<je ~ 4> P 
for Tl p /v n <C 1; AE/cte ~ £ for Tl p /v n 1. For AE/cte > 1 the production process can 
discriminate between different neutrino eigenstates, leading to a loss of their coherence. 

It is interesting to note that for decay of pointlike parent particles at rest (vp — 0) the 
production decoherence parameters vanish: the parameter £ vanishes because it is propor- 
tional to Vp, whereas <j) p is essentially zero because the effective size of the neutrino source 
l p is negligibly small for decay at rest. Physically, the reason for a perfect production co- 
herence in the case v n — > (and good spatial localization of the neutrino detection process) 
is that the oscillation baseline is fixed and, in particular, is independent of the time when 
the parent pions decay. The oscillation baseline is not, however, fixed if the wave packets of 
the parent particles are of finite size; in this case the averaging of the oscillation phase over 
the widths of these wave packets can lead to decoherence effects. 

In the second approach considered in the present paper it is assumed that the position 
of each neutrino emission point is exactly known. Each individual neutrino production 
event is assumed to be fully coherent, and the averaging of the production coordinate over 
the neutrino source is performed simply because the source is extended and the neutrino 
emission can occur in any place inside it. The averaging in this case is done at the probability 
level and is described by the parameter <p p for Tl p /v n <C 1 and by £ for Tl p /v n 3> 1. These 
parameters are just related to the effective size of the neutrino source, which in these limits 
is, respectively, l p and /decay = v n /T. Actually, there is no notion of production decoherence 
in this case, and the suppression of the oscillatory terms in the probabilities is merely due 
to the averaging of the standard oscillation probability over the macroscopic sizes of the 
neutrino source and detector. 

As we have discussed in detail, the two approaches turn out to be equivalent if the 
detection process is perfectly localized in space and time and in addition the spatial size of 
the wave packet of the parent particle can be neglected. This conclusion does not change if 
the particles accompanying the neutrino production (such as muon in tt — > [iv decay) interact 
with the medium or are directly detected. However, if the detection of the accompanying 
particles is used for neutrino tagging, i.e. allows one to establish the coordinate of the 
neutrino emission point in the source with an accuracy a Xfl , this can affect the oscillation 
probabilities by reducing the averaging effects. In the case when the coherent amplitude 
summation does not reduce to the probability summation (e.g., for finite-size pion wave 
packets), this happens due to a better localization of individual neutrino production points 
(assuming that a Xfl < l p , /decay), which improves the production coherence. In the case when 
the use of the probability summation is legitimate, neutrino tagging through the detection 
of the accompanying particles would just mean that the neutrinos are emitted not from the 
whole neutrino source, but from a region of it of the effective size a XfJi . The integration over 
the coordinate of the neutrino production should then be carried out only over this region. 
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In conclusion, we have identified the condition under which coherent and incoherent 
summations over points in the neutrino production region yield different oscillation results. 
This condition can be succinctly stated as follows: the two approaches lead to different 
results whenever the localization properties of the parent particles at neutrino production 
and of the detection process are such that they prevent the precise localization of the point of 
neutrino emission. The difference in the oscillation results is negligible for present accelerator 
experiments, and therefore the standard averaging of the oscillation probabilities over the 
finite spatial extensions of the neutrino source and detector properly takes decoherence 
effects into account. Whether it is possible to devise a realistic experiment that could probe 
this difference remains to be seen. 

The authors are grateful to Thomas Schwetz and Manfred Lindner for useful discussions. 
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